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We use inelastic neutron scattering to explore the evolution of the low energy spin dynamics in 
the electron-doped cuprate Pro.8sLaCeo.i2Cu04_i (PLCCO) as the system is tuned from its non- 
superconducting, as-grown antiferromagnetic (AF) state into an optimally-doped superconductor 
(T c « 24 K) without static AF order. The low temperature, low energy response of the spin exci- 
tations in under-doped samples is coupled to the presence of the AF phase, whereas the low-energy 
magnetic response for samples near optimal T c exhibits spin fluctuations surprisingly insensitive to 
the sample temperature. This evolution of the low energy excitations is consistent with the influence 
of a quantum critical point in the phase diagram of PLCCO associated with the suppression of the 
static AF order. We carried out scaling analysis of the data and discuss the influence of quantum 
critical dynamics in the observed excitation spectrum. 



PACS numbers: 74.72.Jt, 61.12.Ld, 75.25.+Z 

I. INTRODUCTION 

The families of high-transition-temperature (high-T c ) 
superconductors are fundamentally composed of two- 
dimensional copper oxygen planes into which charge car- 
riers, either holes or electrons, are doped. Prior to 
doping charge carriers, the parent compounds of the 
high-T c copper oxides are antiferromagnetically ordered 
insulators whose spin dynamics are well modeled by 
a two-dimensional Heisenberg antiferromagnetic (AF) 
formalism with an anomalously large nearest-neighbor 
exchange coupling (J > 100 mcV) between the Cu 
sites within the CuC>2 pian o 1 ' 2 ' 3 . As spin fluctuations 
may play a crucial role in the mechanism of high-T c 
superconductivitjf 4 * 5 ^, it is imperative to have a compre- 
hensive picture on how the spin dynamics of the undoped 
AF parent compounds evolve as they are tuned towards 
optimally-doped superconductivity. While a comprehen- 
sive picture of this carrier induced modification to the 
spin dynamics has emerged for different classes of hole- 
doped high-T c materials 6 * 7 ^ 8 ^^, experiments exploring 
spin fluctuations in electron-doped copper oxides are just 
beginning 11 . As a consequence, studies of electron-doped 
materials provide a unique litmus for testing the electron- 
hole symmetry in spin dynamical properties. If spin fluc- 
tuations are fundamental to the mechanism of high-T c su- 
perconductivity, they should have universal features for 
all copper-oxide systems. 

In the case of hole-doped high-T c superconductors such 
as YBa2Cu306+^ (YBCO), the most prominent feature 
in its spin excitation spectrum is a sharp magnetic exci- 



tation termed 'resonance' observed by inelastic neutron 
scattering. The resonance is centered at the AF ordering 
wavevector Q = (1/2, 1/2) in the two-dimensional recip- 
rocal space of the Cu02 planes (see inset of Fig. la) 
and is intimately related to superconductivity 5,6,7 . Our 
recent discovery of the resonance mode in the electron- 
doped Pr .88LaCe .i2CuO4_,5 (T c = 24 K) demonstrates 
that the resonance is a universal feature of high-T c cop- 
per oxides regardless of carrier typei^. This unifying fea- 
ture in the spin excitation spectrum of both hole- and 
electron-doped cuprates contrasts the prevalent asymme- 
try known between the commensurate low energy spin 
fluctuations centered at Q = (1/2,1/2) in the electron- 
doped cuprate systemsii* 1 ^ 3 - and the observed incom- 
mensurate spin excitations at Q = (0.5 ± S, 0.5 ± 5) 
in several classes of hole-doped materials 8 ^^ 5 -. While 
the doping dependence of the resonance excitation has 
been observed to follow E r « 5.8ksT c regardless of car- 
rier typei£ and the doping evolution of the incommen- 
surate spin fluctuations in hole-doped materials is also 
controlled by Tji 4 * 1 ^, the doping dependence of the low 
energy commensurate spin fluctuations in electron-doped 
materials remains unexplored. 

For this reason, we chose to investigate the evolu- 
tion of low frequency spin fluctuations in the electron- 
doped cuprate Pr .88LaCeo.i2Cu04_i (PLCCO). Com- 
pared with the prototypical electron-doped copper oxide 
Nd2- a; Ce 2 ;Cu04_5 (NCCO), studying spin fluctuations 
in PLCCO offers several distinct advantages: First, crys- 
talline electric field (CEF) levels of Pr 3+ in the tetrag- 
onal unit cell of PLCCO have a nonmagnetic, singlet 



2 



a 




0.12 0.25 

Ce concentration (a) 

b 




-0.12 J 1 1 1 — 1 1 LI — I 1 

5 10 15 20 25 30 35 

Temperature (K) 



FIG. 1: Phase diagram of Pro.88LaCeo.i2Cu0 4 -,5 (PLCCO) 
and magnetic bulk susceptibility measurements for PLCCO 
samples studied, a) Phase diagram of PLCCO from Kang 
el aim— with the T c = 23 K sample added. The blue ar- 
rows denote the sample compositions and locations on the 
phase diagram that were probed in the neutron experiments. 
The inset shows the direction of Q-scans with the circle de- 
noting (0.5, 0.5, 0) and the cross indicating where the offset 
background was measured—, b) Bulk susceptibility measure- 
ments showing the superconducting phase transitions for the 
samples studied. Both the onset temperature of supercon- 
ductivity and the point of 90% of the normal state response 
(shown as dashed lines) are reported. Susceptibilities have 
been normalized to one another for plotting purposes. 

ground state, which avoids complications arising from the 
magnetic ground state of Nd 3+ in NCCGiS*— Second, 
PLCCO can be tuned from an as-grown nonsupercon- 
ducting (NSC) antiferromagnet to an optimally doped 
superconductor, without static AF order, through an 
annealing processi2ii2i2S; whereas static AF order coex- 
ists with superconductivity in NCCO even at optimal 
doping 1 ^ 21 ! 22 . Finally, the cubic (Pr,La,Ce) 2 3 impurity 
phase arising from the annealing process 2 ^ in PLCCO has 
a nonmagnetic ground state as compared to the mag- 
netic ground state of the (Nd,Ce)2 03 impurity phase in 



NCCO 23 . In this article, we present a comprehensive 
study of the low energy spin dynamics in PLCCO from 
0.5 meV < hcu < 5 meV in a variety of samples as PLCCO 
is tuned from a NSC antiferromagnet into an optimally 
doped superconductor (T c — 24 K) through the anneal- 
ing processings. 

Electron-doped copper oxides differ from the their 
hole-doped counterparts in that they require a post- 
growth annealing treatment in a low-oxygen at- 
mosphere to remove excess oxygen and achieve 
superconductivity— This means that the phase diagram 
of electron-doped copper oxides is three-dimensional as 
a function of both Ce and oxygen concentration (Fig. 
la). There are therefore two distinct ways to traverse the 
phase diagram of PLCCO: samples can be prepared at 
a fixed annealing condition but grown with variable Ce 
doping levels^, or alternatively, samples can be grown 
with a fixed Ce concentration but with variable anneal- 
ing treatmen^i2iS&. For our studies, we utilize the lat- 
ter method with a fixed Ce concentration of x = 0.12 
and the resulting phase diagram is shown in Fig. la. 
It is clear that the annealing process induces a rapid 
suppression of the long-range AF order and the eventual 
emergence of a superconducting phase transition. With 
continued oxygen removal, the superconducting phase is 
enhanced to an optimum T c where the static AF or- 
der is completely suppressed. For samples which ex- 
hibit a superconducting phase, an additional quasi-two- 
dimensional spin density wave (SDW) order appears at 
the disallowed three-dimensional AF ordering wave vec- 
tor Q= (1/2,1/2) in the Cu0 2 planei^. The SDW 
order has an onset of approximately T/v for underdoped 
superconducting PLCCO samples. 

The rest of this article is organized as follows: The next 
section describes the experimental procedure, including 
details of sample preparation and neutron spectrometer 
setup. In section III, we discuss low energy magnetic ex- 
citations in the as-grown, AF ordered NSC PLCCO. Sub- 
sequent sections (IV- VI) cover spin fluctuations in super- 
conducting samples annealed gradually toward optimal 
superconductivity. As AF order is suppressed in the sys- 
tem, the low energy excitations transition from regimes 
coupled to the onset of the AF phase into a virtually tem- 
perature independent regime similar to those observed in 
several heavy Fermion systems known to be near a quan- 
tum critical point (QCP) in the phase diagram 25 i 26 i 27 i 28 . 
In section VII, we discuss the applicability of quantum 
critical scaling in describing the spin dynamics of PLCCO 
with different transition temperatures. Finally, in sec- 
tion VIII we briefly summarize the key conclusions of 
the work. 



II. EXPERIMENTAL 

For our experiments, we grew high quality single crys- 
tal PLCCO samples in an infrared mirror image furnace 
using the traveling solvent floating zone technique. All 
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samples were confirmed to have a mosaic of < 1°. Follow- 
ing their growth, samples were annealed in either a high 
vacuum environment (P < mbar) or in an argon 

gas environment at variable temperatures. The resulting 
magnetic susceptibility and superconducting phase tran- 
sition were measured in a SQUID magnetometer and are 
shown for each sample in Fig. lb. The superconduct- 
ing transition temperature is indicated by both the onset 
of the diamagnetic signal, T c (onset), and through the 
dashed lines showing the 90% position of normal-state 
bulk susceptibility, T c (90%). For the remainder of the 
paper, we will reference these samples by their respective 
onset temperature of superconductivity, T c (onset). The 
as-grown sample, with a mass of —3.0 g, is nonsupercon- 
ducting and is labeled as NSC. The masses of the other 
PLCCO's are -1.8 g, -4.8 g, and -3.0 g for the T c = 21 
K, T c — 23 K, and T c = 24 K samples, respectively. 
An argon annealing atmosphere was used in treating the 
T c = 21 K sample at 940 °C for 24 hoursiS. The other 
two samples were annealed for four days in a high vac- 
uum environment at 765 °C and 775 °C for the T c = 23 
K, and T c — 24 K samples, respectively. 



We reference positions in reciprocal space at wave vec- 
tor Q= {q x ,qy,qz) in A -1 using (H,K,L) (r.l.u.) nota- 
tion, where (H,K,L) = (q x a/2ir,q y b/2w,q z c/2ir) for the 
tetragonal PLCCO unit cell (space group: IA/mmm a = 
b = 3.98 A, c = 12.27 A). Neutron scattering experiments 
were performed at the NIST Center for Neutron Research 
on the SPINS and BT-9 triple-axis spectrometers. Data 
was collected on the cold-neutron spectrometer SPINS 
with a fixed final energy of Ef = 5.0 meV and colli- 
mations of open-80'-sample-80'-open-detector for NSC, 
T c = 21 K, and T c = 23 K samples, while an Ef = 3.7 
meV with identical collimations was used for SPINS ex- 
periments on the T c = 24 K samples. For data collected 
from BT-9, an Ef = 14.7 meV was used with collima- 
tions of 40'-60'-sample-80' -open-detector. Unless oth- 
erwise stated, PLCCO samples were aligned within the 
[H, K, 0] scattering plane for experiments on SPINS, and 
within the [H, H, L] scattering zone for BT-9. Alignment 
in the [H, K, 0] scattering zone allows for an effective in- 
tegration along the c-axis of the quasi- two-dimensional, 
inelastic scattering from the Cu spins in uncorrelated 
CuC>2 planes (through relaxed out-of-plane resolution), 
whereas alignment in the [H, H, L] scattering zone al- 
lows an accurate determination of the onset of the three- 
dimensional AF phase in the system. Masses for the 
samples used in these neutron experiments correspond to 
the masses stated above in the bulk magnetization mea- 
surements with the exception of experiments probing the 
T c = 24 K system. For experiments on this optimally- 
doped PLCCO system, a set of three T c = 24 K samples 
were co-aligned with a total mass of —9 gi£. All samples 
were mounted and loaded into a liquid-He cooled cryo- 
stat, and the experiments were performed in the range 
from T = 2 K to 220 K. 
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FIG. 2: Magnetic order parameters and AF Bragg peaks for 
both the NSC and T c = 23 K samples, a) Order parameter for 
the Neel phase in the NSC PLCCO. All points are integrated 
intensities of the Q=(0.5, 1.5, 0) AF Bragg reflection, b) 
Representative elastic Q-scans through Q=(0.5, 1.5, 0) along 
the [H, 1 + H, 0] direction both below Tm at 40 K and above 
T N at 218 K. c) Magnetic AF order parameter for the T c = 23 
K sample. The sample is aligned in the [H, H, L] scattering 
zone. Points are the peak intensities collected at the Q = 
(0.5, 0.5, 1) Bragg position, d) Elastic Q-scans through the 
(0.5, 0.5, 1) position along the [0.5, 0.5, L] direction. A weak 
AF Bragg reflection appears below Tjv at 2 K and disappears 
above Tjv at 30 K. 



III. AS-GROWN, NSC PLCCO 

We first opted to study the behavior of as-grown 
PLCCO, prior to any annealing treatment. For this 
untreated system, long-range AF order persists up to 
a Tjv = 210 K as shown in Fig. 2a. The rapid in- 
tensity increase below about 80 K arises from the in- 
duced Pr 3+ moment through Cu 2+ -Pr 3+ interaction^. 
Scans through the allowed three-dimensional AF order- 
ing wave vector Q= (0.5,1.5,0) show the appearance 
of a resolution-limited AF Bragg reflection centered at 
Q= (0.5,1.5,0) which diminishes above Tjv (Fig. 2b). 
Due to the noncollinear structure of the zero field spin 
arrangement in PLCCO, the in-plane AF Bragg reflection 
at Q = (0.5,0.5,0) is disallowed, resulting in a gapped 
excitation spectrum at this wave vector^Si. Constant- 
E scans along the [H, H, 0] direction show this low en- 
ergy gap in Figs. 3b-d. In Fig. 3b, Q scans at T = 2 
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FIG. 3: Low energy spin excitations in NSC PLCCO. a) 
Constant-Q scans taken at T — 2 K. The peak signal was col- 
lected at Q = (0.5, 0.5, 0) (black diamonds) and the nonmag- 
netic background was measured at Q = (0.53, 0.53, 0) (black 
triangles). The dashed line is a linear fit to the energy de- 
pendence of the nonmagnetic background. The instrumental 
energy resolution of FHWM= 0.259 meV is determined from 
vanadium scans and is shown as a solid bar under the elastic 
incoherent peak, b-d) Constant-i5 scans along the [H, H, 0] 
direction at Ttw = 0.5, 1.5, and 4.0 meV and collected at T —2 
K. Peaks are fitted by Gaussians centered at Q = (0.5, 0.5, 0) 
and are resolution-limited. 



K through the Q = (0.5,0.5,0) position at hu> = 0.5 
meV show no magnetic scattering, whereas for hu> > 1.5 
meV clear resolution-limited spin-wave peaks are ob- 
served centered at Q = (0.5,0.5,0). These resolution- 
limited spin-wave peaks are well fit by Gaussians and 
give minimum in-plane correlation lengths of t; m in = 
248 ± 40 A at Hlu = 1.5 meV and £ mi „ = 187 ± 21 A at 
%uj = 4.0 me\4S°i. Now plotting in Fig. 3a the constant- 
Q scans at the Q = (0.5,0.5,0) position, this low en- 
ergy gap becomes clearer. Background points taken at 
Q = (0.53, 0.53, 0) are overplotted with raw data taken at 
the peak Q = (0.5,0.5,0) position, showing the presence 
of the low energy gap to be E gap w 1.25 meV. This gap 
energy along with the enhancement in the observed mag- 
netic scattering for hu> > 5 meV is consistent with pre- 
vious studies of the parent compound Pr2Cu02 in which 
the opening of an in-plane anisotropy gap was observed 



in addition to the lower interplane gapi2i2£. 

The temperature dependence of the low energy spin 
waves in this NSC sample is shown in Figs. 4a-c. For en- 
ergies above the gap, the spin-wave excitations increase 
in intensity with increasing temperature (for T < TV) 
following the Bose population factor [n(ui) + 1] = 1/(1 — 
exp^—Twj/ksT)). Removing the nonmagnetic back- 
ground contributions to the scattering, we obtain 5(Q, ui) 
and determine the imaginary part of the dynamic suscep- 
tibility, x"(Q,w), using S(Q,u>) oc [n(w) + l]x"(Q,u;). 
Now plotting x"(Q, lu) in Figs. 4d-f, the T-independence 
of the dynamic susceptibility for E > E gap becomes 
clear. The peaks in x"(Q,o;) a t these energies remain 
T-independent consistent with the excitations simply fol- 
lowing the Bose statistics expected for spin-wave excita- 
tions. The presence of bose-populated spin wave exci- 
tations in this NSC sample is similar to those observed 
in parent compounds La2Cu04 and P^CuGi 1,2 . Above 
TV, however, the low energy spin gap closes and spin 
fluctuations appear in the Tllo = 0.5 meV channel. This 
is most likely due to the emergence of classical critical 
fluctuations arising from the suppression of AF order at 
T = 210 K, or due to a crossover from three-dimensional 
to two-dimensional spin fluctuations as the weak out-of- 
plane Cu exchange coupling breaks down above TV- The 
overall picture of spin dynamics in this NSC sample is 
well described by spin waves arising from the long-range 
AF order in this system, which now provides a baseline 
for studying how these spin excitations evolve as the sys- 
tem is tuned into superconductivity. 



IV. PLCCO, T c =21 K, Tjv=40 K 

We now turn to study an underdoped, superconducting 
PLCCO sample that is annealed to a T c = 21 K with a 
coexisting TV « 40 HM. The resulting low energy fluctu- 
ations around the AF ordering wave vector for this sam- 
ple are plotted in Fig. 5. Raw scattering intensities from 
energy scans at Q = (0.5,0.5,0) from —2.0 < hu> < 5.0 
meV are shown in Fig. 5a for temperatures above and 
below both T c and Tn as solid circles, whereas the non- 
magnetic background collected at Q=(0.53, 0.53, 0) is 
overplotted as solid triangles. One immediate difference 
appearing between the spectra of this underdoped su- 
perconductor and the NSC sample is the appearance of 
a peak in the low-T dynamic susceptibility along with 
the absence of a low energy spin-gap. Figures 5c-e show 
this gapless peak through Q-scans at 7ia;=0.5, 1.5, and 
4 meV, where the T=2 K response displays a clear en- 
hancement at huj=l.b meV. Solid lines show that Gaus- 
sian fits centered at (0.5,0.5,0) on linear backgrounds 
can well describe the observed peaks in Figs. 5 c-e. The 
T = 2 K excitations exhibit widths much broader than 
the resolution-limited spin-wave peaks of the NSC sys- 
tem (Figs. 3c and 3d). Calculating the Fourier trans- 
forms of these Gaussian peaks yields the minimum in- 
plane dynamic correlation lengths of £ m in = 123 ± 21 A 
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FIG. 4: Temperature dependence of the inelastic neutron 
scattering intensity and the dynamic susceptibility in NSC 
PLCCO. a-c) Raw scattering intensities for Q-scans along the 
[H, H, 0] direction for various temperatures at Tiuj — 0.5, 1.5, 
and 4.0 meV. Peaks in scattering are fitted by Gaussians on 
a linear background. A clear gap in scattering at Tiuj = 0.5 
meV persists above 40 K until the three-dimensional order 
breaks down above Tjv. d-e) Measured dynamic susceptibil- 
ity at various temperatures for Tiuj — 0.5, 1.5, and 4.0 meV. 
Dashed lines are Gaussian fits to T—21S K x"(Q,0.5 meV) 
in panel (d), T=2 K x"(Q> 1-5 meV) in panel (e), and T=2 
K x"(Q,4.0 meV) in panel (f). 



at hui = 0.5 meV, £ m j„ = 165 ± 20 A at hui = 1.5 meV, 
and t; m in = 121 ± 27 A at hui — 4.0 meV. Using the 
resolution-limited widths of the NSC PLCCO as refer- 
ence, the true in-plane dynamic spin correlation lengths, 
£, are calculated to be £ = 220 ± 21 A at hui = 1.5 meV 
and £ = 160 ± 47 A at hui = 4.0 meV for the T c = 21 
K PLCCO. The substantially broader Q-widths in Figs. 
5c-e than those of the resolution-limited spin waves in the 
NSC PLCCO suggest that these excitations cannot arise 
from the classical spin-wave scattering from the three- 
dimensional static Neel ordered phase in the sampl er 11 ' 1 ? . 

Upon warming, the low-T peak in the dynamic suscep- 
tibility vanishes and excitations at hui < 3.5 meV popu- 
late upward until the system is warmed above TV. Above 
Tjv, however, there is a crossover in the magnetic re- 
sponse, and the spectral weight for hw < 3.5 meV begins 
to decrease with increasing temperature. This crossover 
is also shown in Fig. 5b in which the intensity of the 
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FIG. 5: Temperature dependence of low energy spin exci- 
tations measured on the T c = 21 K sample, a) Representa- 
tive constant-Q scans at 2 K , 40 K, and 70 K collected at 
Q = (0.5, 0.5, 0) (circles), where the solid black line shows the 
background scattering collected at Q = (0.53,0.53,0) (trian- 
gles). At T = 2 K, a clear enhancement can be seen around 
Tiuj = 2 meV. The solid bracket underneath the incoherent 
elastic line is the measured energy resolution by a vanadium 
standard, b) The temperature dependence of the Tiuj = meV 
and Tiuj = 0.5 meV scattering at Q = (0.5, 0.5, 0). The arrows 
mark the onset of Tjv and T c . c-e) Q-scans of the T c = 21 
K PLCCO along [H,H,0] for huj = 0.5, 1.5, and 4 meV and 
T = 2 K, 40 K, 70 K and 118 K. Center brackets are instru- 
mental resolutions measured by resolution-limited spin-wave 
peaks from the long-range AF ordered as-grown NSC sample. 
The inset in panel (e) shows Q-scans along the c-axis, where 
solid line is the Cu 2+ form factor squared and the dashed line 
marks the background. 



hui = 0.5 meV excitations and the static hui = meV 
SDW moment are both plotted as a function of tem- 
perature. With increasing temperature from 2 K, the 
hui = 0.5 meV fluctuations increase in intensity until the 
breakdown of the static SDW order near Tjv- On further 
warming, these hui = 0.5 meV fluctuations begin to damp 
until disappearing for temperatures above 118 K (Figs. 
5b and 5c). These excitations are also found to be quasi- 
two-dimensional within the Cu02 layers through experi- 
ments which orient the crystal in the [H, H, L] scattering 
plane. This facilitates [0.5, 0.5, L] scans in the inset of 
Fig. 5e, where the L-dependence of the spin excitations 
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FIG. 6: Temperature dependence of x"(Q,^) a t the Q = 
(0.5,0.5,0) position for hui — 0.5, 1.5, and 4.0 meV. a-c) 
x"(Q,w) for T < Tjv (T — 2 K and 40 K). d-f) x"(Q,w) 
for T > Tjv (T = 70 K and 118 K). In all panels, solid lines 
are Gaussian fits to the T = 2 K dynamic susceptibility at 
the corresponding energy. 



at hu> = 0.5, 1.5, 4 meV is rod-like and simply decays fol- 
lowing the Cu 2+ magnetic form factor (solid line). This 
weak correlation perpendicular to the CuC>2 planes in- 
dicates that spin excitations are uncorrelated along the 
c-axis and reflect the highly anisotropic exchange cou- 
pling between Cu sites within the CuC>2 plane (J > 100 
meV) and the much weaker out-of-plane exchange^. 

For huj > 3.5 meV, the magnetic scattering intensity 
at Q = (0.5,0.5,0) is T-independent from T = 2 K to 
120 K as shown in Figs. 5a and 5e. In contrast, the 
spin scattering at hui < 3.5 meV is strongly coupled to 
the appearance of the AF phase. Hence, this regime of 
T-independent S(Q,£j) provides additional evidence of a 
crossover in the spin dynamics that is now instead re- 
flected in the energy scale of the observed spin fluctu- 
ation spectrum. A changeover in the response of the 
system can also be tested by plotting ^"(Q,^) at vari- 
ous energies for temperatures both above and below Tv 
as shown in Fig. 6. The top panels in Figs. 6a-c show 
the measured dynamic susceptibility below Tjy. From 
these, the population of the lowest energy excitations at 
Too = 0.5 meV can be seen to simply follow Bose statis- 
tics with x"(Q,a;) remaining constant. At hui = 1.5 
meV, there exists a slight decrease in the susceptibility 
upon warming to 40 K, most likely due to the vanishing 
peak in the susceptibility (which is present at 2 K, see 
Fig. 5). The susceptibility at hui = 4.0 meV, however, 
recovers this behavior of following the Bose population 
factor with no T-dependence in x"(Q,^) up to T — 40 
K. For temperatures above Tv, the dynamic susceptibil- 



ity decreases sharply with increasing T for all energies 
studied as shown in Figs. 6d-f. This is simply reflective 
of the decrease in magnetic scattering with increasing 
T (for T > Tv) at hui = 0.5 and 1.5 meV. The de- 
crease in x"(Q,cj) at Tiw — 4.0 meV, instead, arises from 
the observed T-independent S(Q, oj) for temperatures up 
to 120 K. This is in sharp constrast to the hui = 4.0 
meV spin-wave excitations in as-grown PLCCO, where 
the intensity of S(Q, ui) is entirely controlled by the Bose 
statistics (Figs. 4c and 4f). An interesting question then 
arises: How do these regimes of energy and temperature 
(hui, T) that couple to either the AF phase or the param- 
agnetic phase (giving rise to this T-independent response 
in S(Q,ui)) evolve with increased doping? We address 
this question below through both the experimental ob- 
servation (Sections V and VI) and data analysis (Section 
VII). 



V. PLCCO, T c =23 K, Tv=25 K 

We now describe results on an PLCCO sample with 
T c = 23 K (Fig. lb). Experiments probing the static 
magnetic ordering of this system show an AF Bragg re- 
flection at the Q = (0.5, 0.5, 1) position with a Tv ~ 25 
K (Fig. 2c). Q-scans through this AF ordering wave 
vector along the [0.5, 0.5, L) direction reveal that this 
weak reflection disappears for T > 25 K (Fig. 2d). In 
Fig. 7a, the measured scattering at Q = (0.5, 0.5, 0) for 
— 1.5 < hui < 5.0 meV is plotted as crossed-boxes, while 
nonmagnetic background collected at Q = (0.56, 0.56, 0) 
is plotted as solid triangles. Examining these low-T fluc- 
tuations, substantial differences appear between the spec- 
tra of this T c = 23 K sample and those of the T c = 21 K 
system. The excitations at the Q = (0.5, 0.5, 0) position 
remain gapless; however, they lack a clearly defined peak 
in the low-T susceptibility at T = 2 K. A comparison of 
the T = 2 K magnetic excitations at hui — 0.5, 1.5, and 
4.0 meV is shown between the T c = 23 K and T c = 21 K 
PLCCO samples in Figs. 7 c-e. The dashed lines show 
Gaussian fits to the scattering observed in the T c = 21 
K sample whose enhancement at hui = 1.5 meV con- 
trasts the continual decrease of scattering intensity with 
increasing energy transfer seen in the T c = 23 K PLCCO 
(crossed-box symbols). Additionally, these plots show 
that the excitations observed in the T c = 23 K system 
have broadened in Q at all energies relative those ob- 
served in the T c = 21 K sample. The Fourier transforms 
of the Gaussian fits in Figs. 7c-c give minimum in-plane 
dynamic spin correlation lengths of £ m i„ = 86 ± 10 A at 
hui = 0.5 meV, £ m j n = 87 ± 12 A at hu> = 1.5 meV, and 
£,min = 67 ± 11 A at hui = 4.0 meV. Correcting for the 
instrumental resolution yields £ = 93 ± 14 A at hui = 1.5 
meV and £ = 72 ± 13 A at hco = 4.0 meV for the T c = 23 
K PLCCO. 

The temperature dependence of the magnetic excita- 
tions in Fig. 7 a shows a similar type of crossover in 
the spin dynamics to those observed in the under-doped 
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FIG. 7: Temperature dependence of the low energy spin 
fluctuations in the T c = 23 K PLCCO. a) Constant-Q 
scans at Q = (0.5,0.5,0) and 2 K, 30 K, 60 K and 100 
K. The solid black line shows the background scattering at 
Q = (0.56, 0.56, 0) and the horizontal bar beneath the inco- 
herent elastic peak is the instrumental resolution, b) Temper- 
ature dependence of the inelastic (hui = 0.5 meV) scattering 
at Q = (0.5, 0.5, 0) and elastic (hu — meV) scattering at 
Q = (0.5, 0.5, 1) Bragg position also shown in Fig. 2c demon- 
strating a magnetic order at T/v ~ 25 K. c-e) Q-scans of the 
T a = 23 K PLCCO along the \H, H, 0] direction for Ulu = 0.5, 
1.5, and 4 meV, respectively, at T = 2 K. The dashed lines 
show identical scans from the T c — 21 K sample with its 
(1,1,0) Bragg intensity normalized to that of the T c — 23 
K sample. Center brackets are instrumental resolutions and 
solid lines are Gaussian fits on flat backgrounds. 



T c = 21 K system. The regimes of these two types of spin 
dynamics are again determined by the energy and tem- 
perature scale of the AF order in the system, now with 
a Tjv ~ 25 K. For T < Tv, there exists a slight enhance- 
ment upon warming of the spin fluctuations at foui < 2.5 
meV, whereas for T > Tv and Uoj > 2.5 meV the ob- 
served scattering intensity remains T-independent over a 
broad temperature range (T = 2 K — > 60 K). The relative 
intensity of the magnetic scattering at T = 100 K is diffi- 
cult to determine due to uncertainties in the nonmagnetic 
background contributions, whereas there was no observed 
change in background scattering between T = 2 K and 
60 K. The enhancement coupled to the AF order is sig- 
nificantly damped and no longer strictly follows the Bose 
population factor. Instead, a weak increase in the popu- 
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FIG. 8: Temperature dependence of the low energy spin fluc- 
tuations and dynamic susceptibility in the T c = 24 K sam- 
ple, a-c) Constant-15 scans at T = 2, 30, 55, 80, and 95 K 
through the Q = (0.5, 0.5, 0) position along the [H, H, 0] di- 
rection. Solid lines are Gaussian fits on linear backgrounds. 
The T = 2, 30, and 55 K data are replotted for compari- 
son from Ref. 11. Solid brackets show the resolution of the 
spectrometer. The increased background scattering with in- 
creasing temperature at %u> = 0.5 and 1.5 meV may arise 
from single/multi-phonon scattering and/or air scattering, d- 
f) x"(Q,w) at T = 2, 30, 55, 80, and 95 K. Dashed lines are 
Gaussian fits to the T — 2 K dynamic susceptibility. 



lation of the %u) = 0.5 meV fluctuations is observed (as 
shown in Fig. 7b) as the system is warmed towards Tv. 
The coupling of these low energy fluctuations to Tv as 
a function of T is also significantly broadened (Fig. 7b), 
in contrast to that of the T c = 21 K PLCCO (Fig. 5b). 
On the other hand, the temperature and energy region 
in which the magnetic scattering intensity is tempera- 
ture independent increases, reflecting the much weaker 
AF phase in this sample along with the reduced energy 
scale of Tv. 



VI. PLCCO, T c = 24 K, Tjv < 600 mK 

Turning now to the final optimally doped T c = 24 K 
sample, previous experiments have shown that there ex- 
ists no static AF order co-existing with superconductivity 
in this sample down to 600 mKi£. The low energy ex- 
citations for this system were first reported in Ref. 11 
and are expanded upon in Fig. 8. Figures 8a-c show the 
raw scattering intensities observed for hw = 0.5, 1.5, and 
3.5 meV at various temperatures. The spin excitations 
in these energies are fitted by Gaussians on linear back- 
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grounds and give minimum in-plane dynamic spin cor- 
relation lengths of £ m j„ = 96 ± 15 A at hui = 0.5 meV, 
£ min = 80± 10 A at Two = 1.5 meV, and £ min = 94±24 A 
at fiuj = 3.5 meV at T = 2 Ki£. These widths are within 
error to those observed in the T c — 23 K sample when 
neglecting the slight change in instrumental resolution on 
changing Ef from 5.0 meV to 3.7 meV. This is reason- 
able since the measured Q-widths are appreciably larger 
than the instrument resolution in both geometries 1 ^. 

Upon warming from T = 2 K to T = 30 K, there is no 
change in the measured magnetic scattering intensity at 
all energies down to Hlo = 0.5 meV. Further increase in 
temperature to T = 55 K renders a slight reduction in 
scattering at hu> — 0.5 meV and 1.5 meV, with no change 
measured in the 3.5 meV excitations. At T = 95 K, the 
peak at hu> = 0.5 meV has completely vanished. There 
exists a strong suppression of the hui = 1.5 meV fluctua- 
tions at T = 80 K. In contrast, there are no changes in the 
Tiuj = 3.5 meV fluctuations up to 80 K. Figures 8d-e show 
the measured x"(Q,w) for the same energies reflecting a 
continued decrease in the susceptibility with increasing 
temperature. The dashed lines show Gaussian fits cen- 
tered at (0.5,0.5,0) for the T = 2 K susceptibility to 
highlight the dramatic decrease in x"(Q,w) at different 
energies (Figs. 8d-f). At the highest measured tempera- 
ture (T = 95 K), the system becomes gapped at hui = 0.5 
meV similar to the two other underdoped PLCCO. The 
absence of any regime in which the dynamic susceptibility 
remains constant, reflective of bosonic excitations simi- 
lar to those observed in the T c = 21 K samples, demon- 
strates a drastic deviation from the two-distinct regimes 
of magnetic response observed in the T c = 21 K PLCCO. 
Therefore, low energy spin excitations in PLCCO evolve 
from coupling to the onset of the AF phase in under- 
doped materials to essentially temperature independent 
from 2 K to 30 K for the optimally doped sample. 

Recently, we have discovered that optimally doped 
PLCCO has a resonance^. Similar to hole-doped 
materials^ 5 -, the resonance in electron-doped PLCCO 
increases in intensity below T c and its energy scales 
with 5.8fcsT c forming a universal plot for all supercon- 
ducting copper oxides irrespective of electron- or hole- 
doping 12 . However, in contrast to hole-doped materi- 
als, magnetic excitations below the resonance in electron- 
doped PLCCO form commensurate and gapless scatter- 
ing (Fig. 8). Therefore, the hour-glass shaped dispersion 
in the magnetic scattering of hole-doped superconduct- 
ing materialsS*2&2ii£ may not be a universal feature of 
all high-T c superconductors. Instead, the resonance itself 
appears to be a universal property of the superconduct- 
ing copper oxides. 



VII. DISCUSSION 

The identification of the continuous suppression of 
Tjv as the optimal superconductivity is approached in 
PLCCO as a function of annealing process^ suggests the 
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FIG. 9: uj/T scaling in the dynamic susceptibility of PLCCO 
with Tc = 21 K and Tjv = 40 K. a) The results from a mini- 
mization of x 2 («) for all x"(u, Q) data with 80 K > T > 30 
K and for data with T < 30 K and Tloj > 2 meV. The bin size 
of 0.5 along the \og(Lu/T) scale is the only assumption made 
in obtaining a = 1. b) The log plot of x"T as a function of 
uj/T shows three distinct regions: 1) data with 80 K > T > 
30 K and with T < 30 K and Huj > 2 meV are within the 
quantum critical scaling regime and collapse onto a universal 
curve. 2) the data at 118 K show a clear high temperature 
departure from this quantum critical scaling behavior. 3) the 
data with T < 30 K and hui < 2 meV show a low temper- 
ature departure from the scaling regime. The breakdown of 
scaling at 120 K could be due to uncertainties in determining 
the backgrounds, c) Summary of the ui/T scaling regime, d) 
x"T as a function of ui/T now plotted with only the data in 
the quantum critical scaling regime shown in panel (c). This 
provides a clearer picture of the universal collapse for data 
within the valid scaling regime. 



possibility of a magnetic QCP, regardless of the precise 
nature of the magnetic structure (homogeneous or in- 
homogeneous). To fully establish the existence of such 
a QCP requires the study of magnetic dynamics, which 
we have shown above using inelastic neutron scattering. 
We focus on two signatures of a QCP. First, a QCP is 
accompanied 31 i 32 i 33 i 34 i 35 i 36 i 37 i 38 by a quantum critical re- 
gion at finite temperatures and finite energies - bounded 
below by a scale which gradually goes to zero as the QCP 
is reached - where the dynamics manifests the excitations 
of the QCP. Second, in this quantum critical regime, the 
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FIG. 10: w/T scaling in the dynamic susceptibility ol PLCCO 
with T c = 23 K, T N = 25 K and T c = 24 K (with no coexisting 
AF order), a) A log plot of x"T as a function of u/T for the 
T c = 23 K sample over-plotted on the universal curve from 
Fig. 9b. By normalizing the (1,1,0) Bragg intensities of 
both samples, we find A — 4.5 which indicates that magnetic 
scattering at Q = (0.5, 0.5, 0) is considerably weaker in the 
T c = 23 K PLCCO. This scaling plot reflects the modified 
scaling regimes in this T c = 23 K sample as discussed in the 
text, b) x" T plotted as a function of u/T for the T c = 24 K 
sample again over-plotted on the universal curve from Fig. 9b. 
c) Scaling relation now plotted for all three superconducting 
samples, T c = 21 K, T c = 23 K, and T c = 24 K, and their 
respective scaling regimes showing a universal collapse onto 
the curve from Fig. 9b. 



dynamics are scale- invariant, a particular form of which 
is an u/T scaling^iS; such an u/T scaling has been 
systematically studied in heavy Fermion metals2£i2§i2£ in 
which the existence of a magnetic QCP is not in doubt^. 
In the following, we will determine how our data can be 
described by the standard QCP theory. 

To perform a scaling analysis for the T c = 21 K 
PLCCO, we note that the u/T scaling can be written 
as x"((5,w)T a = F(Q,u/T), where the scaling expo- 
nent a and the scaling function F(Q,u/T) are deter- 
mined through the best-observed collapse of the data 
onto one universal curved. The quantum critical scal- 
ing is different from the critical scattering from classical 
AF second order phase transition as the former is con- 
trolled by T itself while the characteristic energy scale 



u c in the latter case is determined by reduced temper- 
ature t = \T - T N \/T N {\t\ < l) 40 . For the Heisen- 
berg antiferromagnet with Tv = 40 K, w c should fall 
within the quasi-elastic part of the energy spectra and 
thus not contribute to the observed magnetic scattering 
in Fig. 5. By plotting all data above the 30 K energy 
scale (hu > 3 meV) as a function of u/T and minimizing 
X 2 (Ref. 38), we find a = 1 independent of the func- 
tional form of F(Q,u/T) (Fig. 9a). Using this value for 
the exponent leads to a scaling plot of x"(Q,u)T as a 
function of u/T (Fig. 9b). Such a plot also identifies 
the scaling regime, which is bounded at low energy (2 
meV) and temperatures (Fig. 9c). The data within the 
scaling regime in Fig. 9c clearly show a collapse onto a 
universal curve as plotted in Fig. 9d. We fit the latter 
using F(y) oc y/(l + (y/C) 2 ) and find C = 0.44 ± 0.02. 
This collapse of the data over more than two decades of 
u/T strongly suggests the presence of universal dynam- 
ics. If these universal dynamics indeed reflect a nearby 
QCP, the scaling behavior must break down at low ener- 
gies and temperatures, with the cut-off scale determined 
by the distance to the QCP. Indeed, for small u the u/T 
scaling is observed only at temperatures above about 30 
K (Fig. 9b). The fact that this temperature scale is of 
the order of Tv suggests that the cut-off to scaling is con- 
nected with the development of the AF order. Likewise, 
at the lowest measured temperature, T = 2 K, the u/T 
scaling is seen only at Tiu above about 2 meV (Fig. 9b), 
which, within the error bars, is equal to the low-u; cut-off 
temperature multiplied by the universal constant C . 

Now turning to how this scaling regime evolves in sam- 
ples tuned closer to optimal doping, the same type of 
scaling analysis can be performed for the low-energy mag- 
netic spectra of the T c = 23 K PLCCO. Since our neu- 
tron diffraction measurements at Q = (0.5, 0.5, 1) show a 
T/v ~ 25 K (Fig. 2c), this sample should be much closer 
to the magnetic QCP than the T c = 21 K (T N = 40 K) 
PLCCO. Figure 10a shows the summary of the x"(Q> u)T 
of the T c = 23 K PLCCO over-plotted on the universal 
scaling fit to the T c = 21 K sample (dashed line). A 
scale factor, A, was used to normalize the x"(Q,u;) at 
Q = (0.5,0.5,0) for the T c = 23 K PLCCO to that of 
the T c = 21 K PLCCO, which reflects the fact that the 
Q- width in the former is much broader (Figs. 7c-e). Be- 
sides the overall scale factor, we find that all data with 
Tiu > 1 meV and/or T > 25 K fall on the universal curve, 
consistent with the T c = 23 K PLCCO being closer to the 
QCP point. The u/T scaling regime is seemingly corre- 
lated with the gradual suppression of the AF order as 
optimal superconductivity is approached, indicating the 
presence of a magnetic quantum critical point in electron- 
doped superconductors. However, the dynamic in-plane 
spin correlation length as stated earlier is observed to de- 
crease with increasing T c and decreasing T/v, suggesting 
that this QCP in PLCCO cannot arise from criticality 
towards a three-dimensional long-range AF order at 
K2£. Instead, the data are similar to what happens in 
the Li-doped La2Cu04— and in certain spin-glass QCP 
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heavy Fermion systemsSSiSl. Although there is at present 
no comprehensive theory for the expected properties of 
a metallic spin-glass QCP, the reduced magnetic transi- 
tion temperature and the continuous nature of the tran- 
sition (Fig. 2c) imply that the scaling behavior in the 
T c = 23 K PLCCO should persist down to lower energies 
and temperatures than that of the T c = 21 K sample. 



Continuing this analysis on the magnetic fluctuations 
observed in the T c = 24 K sample, an identical scaling 
plot is shown in Fig. 10b. A constant scale factor was 
again used to normalize the scaling fit from the T c = 21 
K sample to overlay with the data from the T c — 24 K 
sample. As the system is tuned closer to the QCP, there 
should now be larger ranges of energies and tempera- 
tures (huj,T) in which the scaling remains valid. Indeed, 
now all energies and temperatures (for which a suitable 
nonmagnetic background was measured) are shown to 
collapse on the same universal curve determined for the 
T c = 21 K sample (dashed line in Fig. 10b). This encom- 
passes the entire range of energy and temperature probed 
from 0.5 > hu> > 4.0 meV and 2 > T > 80 K with the 
exception of a notable divergence at the fiuo — 0.5 meV 
spin fluctuations measured at T = 2 K in Fig. 10b. 
This divergence from the scaling fit for spin dynamics 
with huj < 1.25 meV may result from the system being 
tuned slightly beyond the QCP. In this case, the low en- 
ergy/temperature spin dynamics crossover into the quan- 
tum disordered regime and no longer obey the scaling 
relation of quatum critical excitations. Since the system 
does not exhibit AF order to at least 600 mK, we can- 
not determine how close the sample is to the Tjv = 
K QCP. Nevertheless, the overall trend for the three su- 
perconducting PLCCO investigated is that of increasing 
regimes of validity for the dynamic ui/T scaling as the 
system is tuned closer to the QCP itself. 



Finally, the valid scaling regimes for all three SC 
PLCCO samples are overplotted in Fig. 10c showing a 
universal collapse onto a common function, with the ex- 
ception of the lowest energy excitations at T — 2 K for 
the T c = 24 K sample as discussed earlier. The universal 
collapse signifying lu/T scaling and the systematic evolu- 
tion of scaling regime as AF order is suppressed with an- 
nealing strongly suggest the presence of a magnetic QCP 
in the electron-doped PLCCO. Recent transport and 
optical measurements on electron-doped Pr2- :E Ce :E Cu04 
(PCCO) have identified singular behaviors, which appear 
compatible with the influence of a QCB^i^. Our results 
suggest that such a QCP may have a magnetic origin, 
possibly due to the presence of a spin-glass QCP. We note 
that previous work has shown the presence of a spin-glass 
QCP in the hole-doped cuprates^, thus suggesting that 
this magnetic QCP might be a common feature in high- 
T c superconductors regardless of doped-carrier type. 



VIII. CONCLUSIONS 

We have systematically measured the doping evolution 
of the low energy spin fluctuations in the electron-doped 
cuprate, PLCCO, as the system is tuned from an as- 
grown NSC antiferromagnet into a phase-pure optimally 
doped superconductor. The as-grown, semiconducting 
PLCCO system exhibits gapped low energy spin fluctu- 
ations consistent with those observed in the parent com- 
pound Pr2Cu04. Tuning the system into an optimally 
doped superconductor changes these low energy spin ex- 
citations dramatically. Instead of a gapped spin wave 
spectrum, the superconducting PLCCO samples exhibit 
gapless low energy spin dynamics which exhibit two dif- 
fering regimes of response in (7ki>,T): the first of these 
coupling to the onset of the AF order in the system and 
the second appearing as spin dynamics whose scatter- 
ing is weakly temperature dependent over large temper- 
ature ranges. The regime that couples to the AF order 
in these samples is seen to evolve from T < 40 K and 
Huj < 2.5 meV in the T c = 21 K sample to T < 25 K and 
hcj< 1.5 meV in the T c = 23 K PLCCO. At optimal dop- 
ing (T c = 24 K), where the static AF phase is suppressed 
to below 600 mK, the observed magnetic scattering for 
T < 30 K or hw > 1.5 meV is T-independent, suggestive 
of the possible influence of quantum critical fluctuations 
in the system. Subsequent scaling analysis of the data 
revealed a collapse of appropriate regions of (hu>, T) onto 
a universal curve for all three superconducting PLCCO 
systems studied, thereby providing microscopic evidence 
for a QCP in the electron doped PLCCO. Additionally, 
the first system discussed, T c = 21 K PLCCO, displays a 
peak in the low-T susceptibility centered at Tq ~ 2.0 meV 
which vanishes upon warming, whereas systems tuned 
closer to optimal doping (T c = 23 K and T c = 24 K) 
display only continuously decreasing magnetic response 
with increasing energy transfer. Our experiments de- 
scribed here have provided a systematic investigation 
into the evolution of the low energy spin dynamics in an 
electron-doped copper oxide, PLCCO, thereby providing 
valuable constraints on microscopic theories seeking to 
model the spin excitations as the materials evolve from 
a long-range ordered antiferromagnet into an optimally 
doped superconductor. 
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